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We have derived the Bethe-Salpeter (BS) equations for the collective-mode spectrum of a mix- 
ture of fermion atoms of two hyperfme states loaded into a moving optical lattice. The collective 
excitation spectrum exhibits rotonlike minimum and the Landau instability takes place when the 
energy of the rotonlike minimum hits zero. The BS approach is compared with the other existing 
methods for calculating the collective-mode spectrum. In particular, it is shown that the spectrum 
obtained by the BS equations in an excellent agreement with corresponding spectrum obtained by 
the perturbation theory, while the Green's function formalism provides slightly different results for 
the collective excitations. 
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I. INTRODUCTION 

Ultracold atoms loaded in optical lattices are ideal sys- 
tems to simulate and study many different problems of 
correlated quantum systems because near the Feshbach 
resonance the atom-atom interaction can be manipulated 
in a controllable way by changing the scattering length a s 
from the Bardeen-Cooper-Schrieffer (BCS) side (negative 
values) to the BEC side (positive values) reaching very 
large values close to resonance. On Bose-Einstein con- 
densation (BEC) side of the resonance the spin-up and 
spin-down atoms can form diatomic molecules, and these 
bosonic molecules can undergo a BEC at low enough 
temperature^— We focus our attention on the BCS side 
where BCS superfluidity is expected analogous to super- 
conductivity. 

In what follows we examine the spectrum of the single- 
particle excitations and the spectrum of the collective 
excitations of a balanced mixture of atomic Fermi gas of 
two hyperfme states with contact interaction loaded into 
an optical lattice. The two hyperfme states are described 
by pseudospins a =t, I- There are M atoms distributed 
along N sites, and the corresponding filling factors / = 
M/N are assumed to be smaller than unity. If the lattice 
potential is sufficiently deep such that the tight-binding 
approximation is valid, the system is well described by 
the single-band attractive Hubbard model: 



H = - J E 



(1) 



Here, the Fermi operator tjr ia (tpi,^) creates (destroys) 
a fermion on the lattice site i with pseudospin a and 
ni,a = V ; i (T V'i,CT is the density operator on site i with a 
position vector r^. \x is the corresponding chemical po- 
tential, and the symbol J2<ij> means sum over nearest- 
neighbor sites. J is the tunneling strength of the atoms 
between nearest-neighbor sites, and U is the on-site in- 
teraction. On the BCS side the interaction parameter 



U is negative (the atomic interaction is attractive). We 
assume h = 1, Boltzmann constant fc# = 1, J = 1 and 
the lattice constant a = 1. 

In the case when the periodic array of microtraps is 
generated by counter propagating laser beams with dif- 
ferent frequencies the optical lattice potential is moving 
with a velocity — v (in the laboratory frame) which mag- 
nitude is proportional to the relative frequency detuning 
of the two laser beams. In a frame fixed with respect 
to the lattice potential the fermion atoms flows with a 
constant quasimomentum p = mv, where m is the mass 
of the loaded atoms. For balance Fermi gases the or- 
der parameter field $_/(«) = —\U\ < ipj^(u)?pj^(u) > 
(or $*(u) = —\U\ < m^(u)i/)\^(u) >) in the mean- 
field approximation varies as <F,- oc Ao exp [2«p.rj]. Here, 
the symbol <> means ensemble average, and A is a 
real quantity which depends on the lattice velocity^ 
In a moving lattice the formation of a BCS superflu- 
idity is possible, but due to the presence of quasimo- 
mentum p the superflow can break down. The stabil- 
ity of balance superfluid Fermi gases loaded into a mov- 
ing optical lattice has been recently studied using the 
second-order time-dependent perturbation theory^, and 
the Green's function formalism^ The superfluid state 
could be destabilized at a critical flow momentum via 
two different mechanisms: depairing (pair-breaking) at 
p pb and Landau instabilities at p cr . The depairing takes 
place when the single fermionic excitations are broken, 
while the Landau instability is related to the rotonlike 
structure of the spectrum of the collective excitations. 
The superfluid state becomes unstable when the energy 
of the rotonlike minimum reaches zero at a given quasi- 
momentum. The numerical solution of the number, gap 
and the collective-mode equations shows that at a zero 
temperature the Landau instability appears before the 
depearing mechanism. - 

Generally speaking, there exist three methods that can 
be used to calculate the spectrum of the collective excita- 
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tions of Hamiltonian ([!} in a stationary (or moving) op- 
tical lattice in generalized random phase approximation 
(GRPA). The first approach uses the Green's function 
method^ - — the second one is based on the Anderson- 
Rickayzen equations jiir— while the third employs the 
perturbation technique^ 

Decades ago, the Green's function approach has been 
used to obtain the collective excitations in the exciton 
problem^ and in the s-wave layered superconductors^ 
In both cases we have deal with a system of inter- 
acting electrons, where two different interactions exist: 
the Coulomb interaction v(r) — e 2 /|r| and a short- 
range attractive interaction g(r) between electrons. In 
the GRPA the collective-mode spectrum can be ob- 
tained by the poles of: (i) the vertex T from the cor- 
responding BS equation for T^iii (ii) the two-particle 
Green's function K by solving the BS equation for the 
BS amplitudes j 10 i 12 ' 13 and (iii) the density and spin re- 
sponse functionsi^ It is worth mentioning that in the 
exciton problem the BS vertex equation^ and the BS 
equation K~ x = K^ 1 - I d - I exc (here is the 

free two-particle propagator, and Id takes into account 
the direct interaction between the electrons, while I exc 
describes their exchange interaction^ ), both provide the 
same collective-mode spectrum, while in the case of su- 
perconductivity, the vertex equation obtained in Ref. (TTj 
in for the collective-mode dispersion (see Eq. (14) in Ref. 
11]) is incorrect because the exchange interaction has 
been neglected (the incorrect equation follows from our 
secular determinant @ assuming a stationary lattice and 
neglecting all elements with 7 and m indices). 

Cote and Griffin^ have obtained the collective ex- 
citation spectrum from the poles of the density and 
spin response functions using the Baym and Kadanoff 
formalisnifA£ They have ignored the long-range Coulomb 
interaction keeping in the direct interaction only the lad- 
der diagrams involving the short-range interaction. The 
exchange interaction involves bubble diagrams with re- 
spect to both the unscreened Coulomb interaction v(r) 
and the short-range attractive interaction g(r). Thus, 
Eq. (2.33) in Ref. [9| corresponds to the BS equation 



K -i = K (o)-i _ Idi and Eq _ (2.32) in Ref. [3 corre- 
sponds to the BS equation K = K + KI exc K. The Cote 
and Griffin Green's function formalism was employed in 
Ref. [6| to study the collective-mode spectrum in moving 
optical lattices by obtaining the poles of the density re- 
sponse function. In what follows we will see that the BS 
approach provides collective-mode spectrum which lies 
slightly below the one obtained in Ref. p if the both 
methods use the same gap and chemical potential. The 
most important advantage of the BS formalism over the 
Greens function approach is that within the BS approach 
we can obtain the poles of density and spin response func- 
tions in a uniform manner, i.e. one secular determinant 
provides not only the poles of the density response func- 
tion, but the poles of the spin response function as well. 

The second method^ that can be used to obtain the 
collective excitation spectrum of Hamiltonian (TT]) starts 
from the Andcrson-Rickayzen equations, which in the 
GRPA were reduced to a set of three coupled equations 
and the collective-mode spectrum is obtained by solving 
the secular equation Det\A\ — 0, where 

(\U\~ l +Il,l J y ,l It, m \ 

A= J 7>! p^+I^ J 7 , m . (2) 

V Il,m J^f^m |^| ~l~ Im,m J 

In a stationary lattice, the quantities I and J are defined 
in Sec. II by Eqs. (ITU1) when p ~ 0. As can be seen, 
the BS secular determinant provides an equation for the 
collective-mode dispersion that depends on all four co- 
herence factors and the associated determinant is 4 x 4, 
while the determinant ([2]) does not depend on the fourth 
coherence factor 7, so that the corresponding determi- 
nant is only 3x3. 

Recently, Ganesh et al£ have used a perturbation ap- 
proach which in the case of a moving optical lattice pro- 
vides the following 3x3 secular determinant (see Eqs. 
(B8) in Ref. @, where there is a negative sign in front 
of the sum in the definition of Xo )'■ 
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Since Det\D\ = det\A p \, where 



of moving optical lattice. 
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(4) 



one can say that the perturbation method is a general- 
ization of the Belkhir and Randeria method to the case 



This paper is organized as follows. In Sec. II, we derive 
the BS equations for the dispersion of the collective exci- 
tations of imbalanced mixture of fermionic atoms loaded 
in a moving optical lattice. In Sec. Ill, we compare our 
approach with the other existing methods, and we discuss 
the stability of balance gases loaded in a moving lattice. 
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It turns out that the BS approach and the second-order 
time-dependent perturbation theory^ provide essentially 
the same spectrum of the collective excitations, while the 
spectrum obtained by the Green's function method- dif- 
fers from the corresponding dispersion obtained by the 
BS approach. But, the Green's function method and the 
BS approach, both predict that as the quasimomentum 
p increases, the rotonlike minimum reaches zero before 
the pair-breaking mechanism takes place. 



II. BETHE-SALPETER EQUATION FOR THE 
SPECTRUM OF THE COLLECTIVE MODES 
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A basic approximation in our approach is that the 
single-particle excitations in the system will be calculated 
in the mean-field approximation, i.e. our single-particle 
Green's function in the momentum and frequency space 

G t,t Gft4 



is the following 2x2 matrix G 



The 



the corresponding matrix elements are defined as follows: 



S«(k) 



iuj m - E + (k: p) 



-£L(k;p) 



where 



A p (k) = 
Bp(k) = 



u 2 p (k) u p {k)v p (k) 
u p (k)v p (k) vl(k) 

Wp(k) -M p (k)wp(k 
-up(k)up(k) wp(k) 



The symbol uj m denotes w m = (27r//3)(m+l/2); /3 = T _1 , 
T is the temperature, and m = 0, ±1, ±2, .... 

Here, 



>(k) 



i[l + x(k;p)/£ P (k)] and 



v p (k) = J \ [1 — x(k; p)/i? p (k)] are the coherent 

factors, and we have introduced the following nota- 
tions x(k;p) = §[£(p + k) + £(k-p)], 77(k ;P ) = 
l[e(p + k)-g( k-p)], £±(k ;p) = r,(k;p) ± £ p (k), 
and E p (k) = v/ x 2 (k; p) + A 2 . The tight-binding form 
of the electron energy is £(k) = 2J(1 — J^cosfci,) — [i. 
As can be seen, the single-particle excitations in the 
mean-field approximation are coherent combinations of 
electronlike E + and holelike fermionic excitations. 

The thermodynamic potential in a mean-field approxi- 
mation can be evaluated as a summation of quasiparticlcs 
with energy E± (k) : 

- 1 In ([1 + exp (-/3£7+(k; p))] [1 + exp {(3E_ (k; p))])}. 

In the case of balance population of fermionic atoms the 
gap and particle number equations can be obtained by 
minimizing the free energy F(A, /) = Sl+/i/ with respect 



FIG. 1: Diagrammatic representation of the Bethe-Salpeter 
amplitude, the single-particle Green's function, the direct and 
exchange interactions. 
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FIG. 2: Diagrammatic representation of the Bethe-Salpeter 
equation for the amplitude ^p^(k, Q). 



to and A: 

f = Itf E [" P (k)/(^ + (k;p) + ^(k)/(^_(k;p)] , 



|C/|^/(S_(k;p))-/(i? + (k;p)) 



1 = ^V 



(5) 



2£ P (k) 



where f(x) = [exp + 1] 1 is the Fermi distribution 
function. 

Next, we derive the BS equations for the spectrum of 
the collective modes in a moving lattice. The collective 
mode energies u> p (Q) and the corresponding BS ampli- 
/* p .t(k,Q)\ 

\Ert4(k Q) 

tude 4'p(k, Q) = vj/^t^ Qj m a movm g optical 

V4+(k'Q)/ 

lattice and spin-polarized systems can be derived in a 
similar manner as in a stationary lattice. Our approach 
is based on the fact that the Hubbard model with on-site 
attractive interaction can be transformed to a model in 
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which the narrow-band free electrons are coupled to a 
boson field due to some spin-dependent mechanismi 12 i 13 
We first use the Hubbard-Stratonovich transformation to 
transform the quartic Hubbard term in (JTJ) to a quadratic 
form. Thus, the quartic problem of interacting electrons 
is transformed to a quadratic problem of non-interacting 
Nambu fermion fields coupled to a Bose field and there- 
fore, we arrive at the problem of the linear response of 
many- fermion systems under a weak bosonic field when a 
transition to a quantum condensed phase is possible. We 
have similar problem in exciton physics when the BEC 
of excitons can be described by applying the powerful 
arsenal of quantum field theory to obtain exact equa- 
tions for the single-particle Green's function, mass op- 
erator and the BS equation for the two-particle Green's 
functioni£&i£i As in the exciton problem, the mass op- 
erator is a sum of two terms: the Fock "Ep and Hartree 



T,h terms. The Hartree term is diagonal with respect 
to the spin indices and it generates the exchange in- 
teraction I exc = 8T,h/SG in the BS equation for the 
spectrum of the collective excitations, while the direct 
interaction Id = ST, p / SG in the kernel of the BS equa- 
tion originates from the Fock term. In Fig. 1 we have 
shown the diagrammatic representations of the leading 
terms of the direct and exchange interactions. In this 
approximation the BS equation for the BS amplitude is 
* = K (°) [h + Iexc] where if (0) is a product of two 
single-particle Green's functions in the mean-field ap- 
proximation. In Fig. 2 we have shown the BS equation 
for \I>p' T (k, Q). The other three equations for *p^(k, Q), 
^^(k, Q) and ^-'-(k, Q) are similar to that one. We 
can write all four BS equations as a single matrix equa- 
tion: 



$ P (k, Q) = -!^S p (k, Q, tu> p ) *p(q> Q) + S M p( k > Q. ^p) E $ p(q> Q) 



(6) 



Here, 



/^ P J - 4 ' t ' t) (k,Q, lWp ), X^' t4 ' t) (k ) Q,t £ j p ) 0\ 
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(t,4,t,4)/ 



p (k,Q,tw P ), ^ P t ' t44) (k,Q, lWp ) 
(TrMYOn, o ^ ^(t,t4,t), 

i j£-(4>4iT>4) th. r\ \ j^(4->t44), 



^ t4 ' t ' t) (k,Q, lWp ), Kg' M ' T1 (k,Q, iu p ) 
(k.Q.tWp), X P J - t44) (k,Q, lWp ) 0/ 



,M = 



( K p UUU) {k,Q,iu: p ), K^'™ ) Qc,Q,iu p )\ 
(t>4,4,4) (k)Q) ^ p)) Xp t,t,t4) (k;QjlWp) 

^.■U,t) (k) q Xp t.t,t,t) (k; q lLUp) 



if, 
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\0 K^' u 'U(k,Q,icj p ), K^' tM {k,Q,tu P ) ) 



and the terms \U\D = K^Id and \U\M = K^I exc represent the contributions due to the direct and exchange 
interactions, respectively. The two-particle propagator Kp (k, Q, ioj p ) is: 

K$*M()l, Q, iuj p -> u + i0 + ) = E G p J ( k + Q; lUJ p + ^m)Gp ( (k; ix J j m )\ tUp ^ u+lQ + = 

/(£ + (k;p))-/(£ + (k+Q;p)) , /(£_(k; p)) - /(£_ (k + Q; p)) 



i0+ 



l -[£ + (k+Q;p)-£_(k;p)] + «0+ 



p p u- [£-(k + Q;p)-£_(k;p)]+zO+ 
kl /(£ + (k;p))-/(£_(k + Q;p)) 



(7) 



p p 



-[£_(k + Q;p)- J E+(k;p)]+*0+ 



a,' 



Here i,j,k,l = {t,|} and cj p = (2tt//3)p;p = 0,±1,±2,... is a Bose frequency, A« = A^(k), = S^'(k) , 
Ip 7 ' = A«(k+Q) and B« = B«(k+Q). To solve the BS equations we introduce a new matrix $ p (k, Q) = f $ p (k, Q), 

and (j ^ and er z are the Pauli matrices. By means of ( ! > P (k, Q) we obtain the following 



where T — , 

<7 a + cr 
equations for collective modes: 



i ^ $ p (k, Q) = - M ^ f [5 p (k, Q, - M p (k, Q, kj,)] f" 1 ! ^ $ p (q, Q). 

k k q 

The condition for existing a non-trivial solution leads to the vanishing of the following secular determinant: 
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where the following symbols have been introduced: 



1 



K, 



a.b 
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/ (E_ (k; p)) - / (£+ (k + Q; p)) / (E_ (k + Q; p)) - / (E + (k; p)) 



2N 
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E<Q^ 
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w + fi p (k, Q) - £ P (k, Q) u + fi p (k, Q) + £ P (k, Q) 

/ (E_ (k + Q; p)) - / (£?_ (k; p)) / (£ + (k + Q; p)) - / (E + (k; p))" 



Q 



w + ftp(k,Q)+e p (k,Q) w + fi p (k,Q) -e p (k,Q) 

/(£7_(k+Q;p))-/(£;_(k;p)) , / (£+(k+ Q; p)) - / (E+(k; p)) 



w + ft p (k,Q) +e p (k,Q) 



w + fi p (k,Q)-e p (k,Q) 



Here, e p (k,Q) = £ p (k + Q) + £ p (k), e p (k,Q) = £ P (k + 
Q) - S p (k), fi p (k, Q) = 77 P (k) - ry p (k + Q), and a and 
b are one of the following form factors: 

7k, Q = M k K( k + Q) + v p (k)v p (k + Q), 
l KQ = w P (k)Mp(k + Q) - v p (k)v p (k + Q), 
7 p k, Q = itp(k)w p (k + Q) - v p (k)w p (k + Q), 
m k,Q = M P (k)z;p(k + Q) + Wp(k)u p (k + Q). 

As Q — > in accordance with the well-known Goldstone 
theorem, there exists a solution co — > 0. In this case all 
J, K and L vanishes, and the secular equation reduces 
to the gap equation written as = 1 + |{/|/ 7 =i l7= i. 

At a zero temperature, before the pair breaking sets 
in, we have E + (k;p) > 0, £L(k;p) = — J5 + (k;p), and 
the secular determinant assumes the form 
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(9) 

Here we have introduced symbols and J p & : 



/P = J_ <Q & k,Q £ p( k ^Q) 

JP = i v < Q frE, Q [" + "p(k,Q)] 

Q ' b A^ [ w + fip(k,Q)] 2 - £p (k,Q)' 

and the quantities a p q and b p q could be one of the four 
form factors: ; p q , m p Q , 7 P Q or7 P Q . 

The vanishing of the secular determinant (|9|) provides 
the spectrum of the collective excitations w(Q). I p b 
and J p b will display singularities if w(Q) = £+(ko; p) + 
£+(Q — k ;p) at a particular k . This means that the 
supcrfiuid state is not stable because the Cooper pairs 
break into two fermionic excitations.— For fixed p and 
different Q the expression _E + (k;p) + E + (Q — k;p) is 
bounded from below by the threshold line uj p h (Q) = 
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FIG. 3: Collective-mode spectrum Q(q) in ID optical lattices 
for quasimomentum p — 0.21 calculated by using the secular 
determinant |[9j). We set / = 0.5 and \U\ = 2J. The superfluid 
gap and chemical potential are A = 0.419J and /i = 0.655J. 
The insert shows the results from FIG. 4b in Ref. but 
calculated with A = 0.420J and fi = 0.624J. 

(mink [E + (k; p) + E+(Q — k; p)]. Our numerical calcu- 
lations in ID, 2D and 3D show that the threshold line 
w p l (Q) is above the spectrum of the collective excitations 
w(Q), and therefore, the rotonlike minimum will reach 
zero before the pair breaking sets in, i.e. the superfluid 
state is destabilized due to the Landau mechanism. 



III. COMPARISON WITH OTHER 
APPROXIMATIONS 

Our numerical calculations show that there is an ex- 
cellent agreement between the dispersions obtained by 
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FIG. 4: Collective-mode spectrum Q(q) in ID optical lattices 
for quasimomentum p = 0.12 calculated by using determinant 
©. We set / = 0.9 and \U\ = 23. The superfluid gap and 
chemical potential are A = 0.351J and fj, = 1.697J. The in- 
sert shows the results from FIG. 6b in Ref. [(|, but calculated 
with A = 0.350J and ft = 1.690J. 
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FIG. 5: Collective-mode spectrum Q(q) in 2D optical lattices 
for quasimomentum p = (p,p), where p = 0.6502/7r\/2 calcu- 
lated by using the secular determinant ©. We set / = 0.5 
and \ U\ — 12 J. The superfluid gap and chemical potential are 
A = 4.9699J and fj, = 0.7455J. The insert shows the results 
from FIG. 9d in Ref. [6| with the same A and fi. 

Ganesh et al. in Ref. [j| by using determinant ([3]), and 
by our secular determinant (j9]) not only in the case of 
a stationary lattice (see the conclusion section in Ref. 
[IH), but in the case of a moving lattice as well. This is 
due to the fact that the terms with index 7 in the secular 
determinant © are very small. This is due to the fact 
that the terms with index 7 in the secular determinant 
(|9]) are very small. 

In Fig. 3 and Fig. 4 we have compared the collective 
mode spectrum in ID lattices obtained by the BS ap- 
proach and by the Green's function method^ It is worth 



FIG. 6: Critical p cr (lower curves) and pair-breaking p p b (up- 
per curves) quasimomenta in units [it /a] as a functions of U / J 
in ID lattices and / = 0.9. The dashed curves represent our 
calculations, while the dotted curves are the results presented 
in Fig. 7 in Ref. @|). 

mentioning that in the ID case our chemical potentials 
and gaps obtained from Eqs. are not the same as in 
Ref. [6j, and therefore, the direct comparison between 
the corresponding dispersion curves is not correct, but 
the curves arc quite similar. When using the same values 
for the chemical potential and the gap, our dispersion 
curve is slightly below the corresponding curve obtained 
in Ref. @ (see Fig. 5). 

In Fig. 6 we have presented the critical p cr (lower 
curves) and pair-breaking ppi, (upper curves) quasimo- 
menta in ID lattices as functions of U/J. The Green's 
function and BS methods provide similar results for 
-2.5 < U/J < -1.8, but in the interval -1.7 < U/J < 
—1.1, both the pair- breaking and the critical quasimo- 
mentum curves differ from the corresponding curves pre- 
sented in Fig. 7, Ref. Q. The differences between the 
pair-breaking curves are probably due to the different 
solutions of gap and number equations, but the differ- 
ences between the critical momentum curves are most 
like due to the different techniques employed to obtain 
the collective-mode spectrum. It can be seen that ac- 
cording to our curves the Landau instability is still the 
mechanism which destabilized the superfluid states. 

It is interesting to compare the analytical results for 
the speed of sound in a weak-coupling regime in ID op- 
tical lattices obtained by both approaches. In a sta- 
tionary lattice we follow the calculations by Belkhir and 
Randeriai^ to reduce our secular determinant © to the 
following 2x2 determinant: 

Q 2 q( M 2 -4J M +A 2 ) _ u£a m 



where a = \U\/ (ir^/ 4Jfi — fi 2 . The solution io{Q) = cQ 
provides the speed of sound which in a stationary lat- 
tice is c s = v 7 ! — o^y^Vp — A 2 , where the Fermi veloc- 
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ity is Vf — y/4JfJ. — /i 2 . Using the same steps as in a 
stationary lattice, we have obtained the following 2x2 
determinant in ID moving optical lattices: 



where a = |[/|/(7r^4J 2 cos 2 (p) - (2J - /n) 2 ) and 



an 

uq I {2J-n)Qa tan(p) 
"2A "T" 2A 



" 2A 



{2J—n)Qa tan(p) 
2A 



a 



= 0. 
(12) 



an 



2 a(/x 2 — 4J/i + A 2 ) w 2 a (2 J - y)Qau tan(p) 1 2 2 [2 J 2 - A + 4 J 2 cos(2p)]Q 2 a tan 2 (p) 



4A 2 cos 2 (p) 



4A 2 



2A 2 



Q 2 tan 2 (p) 



2A 2 



(13) 



The last determinant provides the following expression for the phononlike dispersion in the long- wavelength limit: 



c s = (2 J - fi) tan(p) + Vl^aJ Vl - A 2 



1 + - — - tan 2 (p) 



(14) 



V F = sj 4J 2 cos 2 (p) - (2 J - ) u) 2 + 4J 2 [2 - 3cos 2 (p)] tan 2 (p). 
The speed of sound obtained in Ref. |6j is 



c s = (2 J - n) tan(p) + \/l - a^j Uj - A 2 
V F = v/4J 2 cos 2 (p) - (2J-/i) 2 . 



1 — Ct T , , 

1 + — 2— tan 2 (p) 



(15) 



r 



It should be mentioned that in the case when A/ J < 
1, the term 4 J 2 [2 — 3cos 2 (p)] tan 2 (p) is more important 
than the term A 2 (l — a) tan 2 (p)/a in (fT4|) and the term 
A 2 (l - a) tan 2 (p) /a 2 in (JI5J|. 



IV. SUMMARY 

In conclusion, we have derived the BS equation in 
GRPA for the collective-mode spectrum of supernuid 



Fermi gases in moving optical lattices assuming that the 
system is described by the attractive Hubbard model. We 
have compared collective excitation spectrum obtained 
by the BS approach with the corresponding spectrum 
obtained by applying the Greens function formalism and 
by the perturbation theory. The BS results are in an ex- 
cellent agreement with the results obtained by perturba- 
tion theory, while the Greens function formalism provides 
slightly different results. 
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